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We observe the suppression of the ID transport of an interacting elongated Bose- 
Einstein condensate in a random potential with a standard deviation small com- 
pared to the typical energy per atom, dominated by the interaction energy. Numer- 
ical solutions of the Gross-Pitaevskii equation reproduce well our observations. We 
propose a scenario for disorder-induced trapping of the condensate in agreement 
with our observations^. 



1. Introduction 

(~| ' Coherent transport of waves in disordered systems is a topic of primary 

O ■ importance in condensed-matter physics, for example in the description of 

normal metallic conduction, superconductivity and superfluid flow, and has 
relevance also to optics and acoustics^. The presence of disorder can lead 
to intriguing and non- intuitive phenomena such as Anderson localization^, 
^ \ percolation dynamics^, and disorder-driven quantum phase transitions to 

Bose-glass^ or spin-glass® phases. The main difhculty in understanding 
quantum transport arises from the subtle interplay of scattering on the 
potential landscape, interferences and interparticle interactions. Due to a 
high degree of control and measurement possibilities, dilute atomic Bose- 
Einstein condensates in optical potentials have proved an ideal system in 
which to revisit traditional condensed matter problems, and recent theoreti- 
cal works discuss disorder-induced phenomena in this context^. Apart from 
the (undesired) effect of a rough potential on trapped cold atoms on atom 
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chips®, there have been few experiments on BEC in random potentials^' 

In this experiment^, we study the axial expansion of an elongated BEC 
in a magnetic guide, in the presence of a disordered potential generated by 
laser speckle. We observe that a weak disordered potential inhibits both 
the expansion and the centre-of-mass (COM) motion of the condensate. 

2. Laser Speckle — a disordered potential for atoms 

Laser speckle^^ is the random intensity pattern produced when coherent 
laser light is scattered from a rough surface. Such patterns arise from 
interferences between wavefronts coming from different scattering sites. In 
the Fraunhofer limit, a speckle pattern does not depend on the details of the 
scattering surface, but follows a well-defined statistical distribution. The 
intensity distribution is exponential P{I) = cxp (—I/aj) with the intensity 
standard deviation equal to the local mean intensity: c/ = (/(z)). The 
typical distance Az between speckle 'grains' can be characterised by the 
half-width of the autocorrelation function. For a circular aperture, this 
is an Airy function, with the first zero located at Az — 1.22X1/ D, where 
A is the laser wavelength, I the distance from the scattering plate to the 
observation plane, and D the aperture diameter at the scattering surface. 




Figure 1. Left: Optical setup used to create the random speckle potential. The con- 
densate is at the focal point of the lens system, with its long axis oriented perpendicular 
to the page. Right: Example of speckle intensity profile, with condensate to scale. The 
speckle potential is effectively ID for the trapped condensate. 

In our experiment, a blue-detuned laser beam is shone onto the atoms 
through a scattering plate, as shown in Fig. 1. The beam is derived from 
a tapered amplifier, injected by a free-running diode laser at A '--^ 780 nm 
and fibre-coupled to the experiment. The out-coupled beam is expanded 
and then focused onto the condensate, the fibre out-coupler and lenses be- 
ing mounted on a single small optical bench, aligned perpendicular to the 
long axis of the BEC. Inserting the scattering plate in the position shown 
projects an optical speckle potential onto the atoms, with an intensity dis- 
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tribution /(r) which is simply the Fourier transform of the phase distribu- 
tion at the scattering plate. The scattered beam has a total power of up 
to 150 mW and diverges to an rms radius of 1.83 mm at the condensate. 
The mean intensity (the gaussian envelope) of the beam can be assumed 
constant over the region where the atoms arc trapped. 

To calibrate the speckle pattern, the optical set-up is removed from the 
BEC apparatus and the intensity distribution observed on a CCD camera 
at the same distance I as the atoms. Taking images with various beam aper- 
tures D, we verify the exponential intensity distribution P{I), and compute 
the autocorrelation function to obtain the grain size Az. Taking into ac- 
count the modulation transfer function of the camera^^, we find that the 
measured grain size follows the prediction. For our setup, / = 140(5) mm 
and D = 25.4(1) mm, giving Az = 5.2(2) /xm. In this experiment, we pro- 
duce condensates with an aspect ratio of 100, typical Thomas- Fermi half- 
length Ltf = 150 /um and radius i?TF = 1-5 /xm. The trapped BEC occupies 
about 45-50 minima along its length, but experiences an almost constant 
potential in the radial directions. (Along the axis of the laser beam, the 
typical length scale of the speckle grains is much longer, ~ Az^/X = 35 /um.) 
The speckle potential is therefore effectively one-dimensional (ID) for the 
atoms: Rtf < Az <C Ltf- This is important as systems in ID potentials 
are known to be more sensitive to the effects of disorder 

We calibrate the speckle intensity {I{z)) as a function of the fibre- 
coupled laser power, from which the light-shift potential V{z) oc I(z)/S 
is determined^. The longitudinal positioning of the optical set-up at a 
distance I from the atoms is controlled to within 5 mm, leading to an un- 
certainty of 3.5% on the speckle intensity. Since the light is blue-detuned 
from resonance {\S\ > 0.15 nm). the potential is repulsive for the atoms. For 
the laser intensities used in this experiment, the standard deviation of the 
optical potential ay = {V{z)) is always smaller than the chemical potential 
fj, of the condensate. 

Other recent experiments have also generated disordered optical poten- 
tials for BEC: using red-detuned laser speckle^ and by imaging randomly 
patterned structures onto an optical lattice^^. 

3. Experimental Sequence 

Condensate preparation. We generate an elongated (quasi-lD) BEC 

of *^Rb atoms in the 5Si/2|^ = l,Tnp = —1) state, using an iron- 
core electromagnet loffe-Pritchard trap^^ with final trap frequencies of 
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a;_L = 27r X 660(4) Hz radially and = 27r x 6.70(7) Hz axially. We 
obtain condensates of typically 3.5 x 10"'' atoms and chemical potential 
IJ,/2Trh ~ 5 kHz. In such an elongated trapping geometry, phase fluctuations 
can be important^^. However, at our estimated temperature of 150 nK, the 
phase coherence length is of the order of the condensate length. 

The laser speckle is turned on at the end of the rf evaporation ramp 
and we wait a further 200 ms in the presence of the rf field to ensure that 
the condensate is in equilibrium at the end of the sequence. 

Axial expansion. We next open the magnetic trap in the axial (dipole) 
direction while keeping the transverse confinement and the random poten- 
tial unchanged, thus converting the trap into a long, uniform, magnetic 
guide. The condensate expands along this guide, the axial expansion being 
driven by the repulsive interatomic interactions. 

Due to hysteresis in the electromagnet and coupling between the dipole 
and quadrupole poles, remanent magnetizations can produce an axial trap- 
ping potential even when the dipole excitation current is reduced to zero. 
For this reason it is necessary to invert the current in the dipole excita- 
tion coils. To avoid spin-flip losses to non-trapped hyperfine states, the 
minimum magnetic field must not cross zero, thus setting a lower limit for 
the dipole current. By extrapolating the variation of the condensate dipole 
mode frequency with the current, we estimate the final axial trapping fre- 
quency to be less than a;^ ^ 27r x 0.4 Hz, which is compatible with the 
linear expansion observed (sec Fig. 2). The trap opening is carried out over 
30 ms, in order to minimize trap loss and heating. Numerical simulations of 
the 3D Gross-Pitaevskii (GP) equation confirm that this 'slow' relaxation 
of the axial confinement reduces only slightly (by ~ 4 %) the asymptotic 
expansion velocity of the BEC with respect to a 'sudden' trap opening. 

Imaging. After an axial expansion time r, the magnetic fields and laser 
speckle are switched off and the condensate is imaged by absorption after a 
further 15 ms timc-of- flight. From axial profiles of these images, the COM 
position and rms half-length L of the atomic cloud are determined. 

4. Inhibition of Transport 

In Fig. 2a we plot the half-length L and COM position of the condensate as 
a function of axial expansion time r. In the absence of the speckle poten- 
tial, the condensate expands axially, reaching a linear expansion velocity of 
2.47(3) mms~^, consistent with the numerical predictions of the GP equa- 
tion. When a small speckle potential ay = O.27 is applied, the expansion is 
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decelerated and the BEC eventually stops expanding. This effect becomes 
more pronounced at higher speckle potentials, as shown in Fig. 2b. The 
BEC also acquires a constant COM velocity of 5.1(2) mms^^ due to a mag- 
netic 'kick' during the trap-opening. This motion too is decelerated and 
finally stopped in the presence of a weak speckle potential. 

By displacing the speckle potential by small distances in the vertical 
(radial) direction, we can project a completely different realization of the 
speckle potential onto the atoms. Since the BEC in the magnetic trap al- 
ready covers many peaks of the speckle potential, we find that the expansion 
dynamics are fairly insensitive to the exact random potential. 
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Figure 2. a) Rms half-length L and COM position (inset) of the BEC as a function 
of axial expansion time t, for three different amplitudes of the disordered potential: 
cry = 0, 0.2/i, 0.7fi. The axial confinement is relaxed during the first 30 ms (grey band), 
b) Variation of rms half-length L with speckle potential ay for t = 115 ms. The three 
points <C> correspond to the data shown in a). 



5. 'Disorder-induced trapping' 

We are able to interpret the observed behaviour by studying the density 
profiles obtained by numerical integration of the GP equation. By assum- 
ing tight radial confinement, that is hoj± ^ hui^, iJ,,kBT, this equation is 
simplified to ID. The static speckle potential is generated numerically, with 
Az = 0.049Ltf- In Fig- 3 we see that this model reproduces well the trends 
observed experimentally, despite the approximation to ID. 

In the harmonic trap, the BEC is in the interaction-dominated, Thomas- 
Fermi regime: n{z,t = 0) = max(0,/i — O.Smw^z^ — V{z))/gii:). Since the 
BEC healing length ^ is much smaller than Az, the density profile simply 
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follows the modulations of the speckle potential. These modulations are 
not sufficiently strong to fragment the BEC, as cry < /i. After switching 
off the axial confinement, the BEC expands (for ay — 0) according to the 
scaling theory For times lu^t > 1 the fastest atoms populate the wings of 
the BEC, where the density is lowest; conversely, at the centre of the BEC, 
the density and therefore the interaction energy remain relatively high, and 
the kinetic energy very low. In the presence of disorder we can therefore 
expect very different behaviour in each of these two regions. 

Figure 3 shows a simulated density profile obtained for cry = 0.2/i and 
LUzT = 10, at which time the BEC has stopped expanding. Looking more 
closely at the central region — Ltf/2 < z < Ltf/2 (Fig. 3a) we find that the 
condensate density modulations have mainly the same length-scale as the 
speckle potential, and that these modulations are stationary. Since in this 




Figure 3. Top Time evolution of the rms size L of the BEC in the random potential, 
for various speckle amplitudes ay, as obtained by numerical integration of the ID GP 
equation. Bottom Simulated density profile (solid) at time lu^t = 10 during axial 
expansion, with o-y = 0.2^ (grey line represents V{z)/gi£)). Right Zoom on density 
profile (solid) — note the different length scales. The approximation n{z) = max(0, Ji{t) — 
V (z)) / gij:) (dashed) reproduces the density only in the central region (a) (see text). 
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region the interaction energy remains dominant, we can suppose that the 
Thomas-Fermi approximation still holds. By defining an effective chemical 
potential Jl{t), which decreases with the density during expansion, we can 
write the density profile as: n{z,t) = max(0,/l(f) — V{z))/giD- This is a 
good approximation close to z = (see Fig. 3a), showing that the local 
density adiabatically adapts to the speckle potential as 'p{t) decreases. In 
this central region, the expansion stops when the condensate encounters 
two large peaks of the speckle potential with amplitudes larger than 
and the condensate then fragments. 

In the low density regions, for example — S.ILtf < z < — 2.9Ltf 
(Fig. 3c), the density profile is not stationary, and furthermore the length- 
scale of the density modulations is much smaller than Az. By using Fourier 
analysis, wc identify the characteristic length-scale as the healing length of 
the trapped condensate, ^ = h/y/2rnjl. This can be understood by conser- 
vation of energy, implying that the kinetic energy per atom is of the order of 
the typical energy per atom /i, of the initial BEC. In this region, the atoms 
are therefore almost free particles, interacting with the disordered poten- 
tial. The BEC wavefunction undergoes disorder-induced multiple transmis- 
sions and reflections and is ultimately blocked by a high peak V{z) > n 
of the speckle potential. The continually evolving contributions of the low- 
density wings of the BEC are responsible for the fluctuations of L observed 
in the simulation for ay ^ 0.2/i, even once the core of the wavefunction is 
localized. At intermediate distances (Fig. 3b), the Thomas-Fermi approx- 
imation is no longer valid, and the density profile exhibits time-dependent 
modulations with a length scale intermediate of ^ and Az. 

To understand the role of the disordered potential in this trapping sce- 
nario, it is useful to compare this situation with that of a periodic potential, 
with a lattice spacing Az and depth Vq = 2cry. Our model predicts similar 
behaviour in the central region, but differs in the low density wings, where 
the condensate would continue to expand due to tunnelling between the lat- 
tice sites. In the lattice, the condensate fragments when the central density 
reaches the value no — Vo/gm, independent of the lattice spacing. In the 
case of the disordered potential, this final density depends on the statistical 
distribution of the optical potential. We can calculate the probability of 
finding two speckle peaks of a given height within a given distance, which 
when combined with the condensate expansion dynamics leads to the follow- 
ing estimate for the final peak density: no — 1.25(cry/(7i£)) ln(0.47LTF/Az). 
This formula, dependent on both ay and Az, is in good agreement with our 
numerical findings and will be the subject of future experimental work^^. 
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6. Conclusion 

We have investigated the transport properties of a BEC in a disordered 
potential, observing a cross-over from free expansion to absence of diffusion 
as the strength of the disorder increases. Our experimental findings are 
supported by numerical simulations and we have discussed a theoretical 
model which describes a scenario for disorder-induced trapping. In the 
future it would be interesting to further investigate this highly controllable 
system, for example by changing the correlation length of the disorder or 
using Bragg spectroscopy to probe the momentum spectrum of the BEC. 

This work is supported by the CNRS, Ministcrc dc la Recherche, DGA, 
EU (IST-2001-38863, MRTN-CT-2003-505032), ESF (BEC2000+), Marie 
Curie Fellowship (J.R.), IXSEA-OCEANO (M.H.) and INTAS (211-855). 
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